РСС/ОСЕ RS79-01 


INVESTIGATION OF METHODS 
FOR CONVERTING THE FCC GROUND WAVE 
FIELD INTENSITY CURVES TO 
THE METRIC SYSTEM 


BY JOHN H. McMAHON 

RESEARCH & STANDARDS 
DIVISION 

OFFICE OF CHIEF ENGINEER 


JANUARY 1979 


FEDERAL COMMUNICATIONS COMMISSION 
WASHINGTON, D.C. 20554 


FCC/OCE RS 79-01 


INVESTIGATION OF METHODS FOR CONVERTING THE FCC GROUND 
WAVE FIELD INTENSITY CURVES TO THE METRIC SYSTEM 


RESEARCH AND STANDARDS DIVISION 
OFFICE OF CHIEF ENGINEER 


JANUARY 1979 


FEDERAL COMMUNICATIONS COMMISSION 
WASHINGTON, D.C. 20554 


Federal Communications Commission 
Office of Chief Engineer 
Research and Standards Division 


8579-01 


INVESTIGATION OF METHODS FOR CONVERTING THE ЕСС GROUND 
WAVE FIELD INTENSITY CURVES TO THE METRIC SYSTEM 


by 
John H. McMahon 


Washington, D.C. 20554 


January 30, 1979 


As 


TABLE OF CONTENTS 


SUBJECT 
List of Figures 
Abstract 
Introduction 
Survey of Methods for Predicting Ground Wave Field Intensities 


Computer Programs Developed for Ground Wave Field Intensity 
Predictions. 


Description of Computer Program Listing 
Discussion and Recommendations 
Acknowledgment 

Bibliography 

Figure 1 

Figure 2 

Figure 3 

Figure 4 


Listing of Fortran Programs used in Report for Calculation of 
Field Intensities 


РАСЕ; 


15 


16 


17 


18 


-2- 


LIST OF FIGURES 


TITLE 


Figure 1, Comparison of FCC and Computer Predicted Ground Wave 
Field Intensities, .55 MHz, Ground Conductivity 2 millimhos, 
Relative Dielectric constant 15. 


Figure 2, Comparison of FCC and Computer Predicted Ground Wave 
Field Intensities, .55 MHz, Ground Conductivity 5000 millimhos, 
Relative Dielectric constant 80. 


Figure 3, Comparison of FCC and Computer Predicted Ground Wave 
Field Intensities, 1.6 MHz, Ground Conductivity 2 millimhos, 
Relative Dielectric constant 15. 


Figure 4, Comparison of FCC and Computer Predicted Ground Wave 
Field Intensities, 1.6 MHz, Ground Conductivity 5000 millimhos, 
Relative Dielectric constant 80. 


Abstract 


This report describes various methods given in the literature for 
calculation of ground wave field intensities in the Standard Broadcast 

band. ‘This examination was made for the purpose of developing a computer 
program which could be used to revise the present FCC Ground Field Intensity 
Curves to the metric system. No method of calculation was found which 
exactly duplicated field intensities given in the existing FCC curves. 

The report contains four figures comparing field intensities predicted 

by several methods and by existing FCC curves. Also included in the report 
is a listing of the computer programs developed during the study for calcula- 
tion of field intensities. 


Йо 


INTRODUCTION 


It appears desirable in keeping with United States Government directives 
to replace English units and, accordingly, to revise the FCC Ground Wave 
Propagation Curves contained in Part 73 of the Commission's Rules to show 
distances in kilometers rather than miles. This conversion appears trivial 
since it would seem only necessary to divide the mileages by 1.609 and 
change the distance scale to kilometers. However, to be useful, major 
graphical divisions should be assigned to integral numbers of kilometers 
to simplify interpolation by the users. Accordingly, this revision requires 
at a minimum, a re-drafting of these curves. If it is desired only to 
revise the distance scales, an inverse field of 160.9 Millivolts per 

meter should be used. If it is deemed necessary to use 100 Millivolts 

per meter at 1 kilometer, complete re-drafting of these curves is required. 


The FCC Broadcast Bureau has pointed out that the present Ground Wave curves 
contain some original drafting errors which should be corrected during the 
metric revision. Since it also appears desirable to develop a computer рго- 
gram to calculate ground-wave field intensities in the Standard Broadcast 
band, the logical method to be used in the revision of these curves is to 
first develop the computer program. This program could be used then with 
available computer graphic programs to redraw the curves. This computer 
program would also be available for direct calculation of ground-wave field 
intensities. 


An investigation has been made of the feasibility of developing a computer 
program for this purpose. Because it is not desired to perpetuate any 
errors in the present curves, two simple methods of developing such a 
program cannot be used. These methods are: 1. Develop a computer program 
which uses polynomial expressions fitted to the present curves and which 
contains the necessary logic to select the proper expression for the given 
frequency and ground constants. 2. Digitize the curves into sets of tables 
and devise table-look programs and interpolation programs to make the proper 
table selections for frequencies and ground constants. With these two 
methods disqualified, the programming method selected should use basic 
mathematical expressions relating field intensities to frequency, ground 
constants, and effective earth curvature. 


A survey of the literature has been made, and three different computer 
programs incorporating different methods of calculating ground wave field 
intensities have been developed. Fields predicted by these programs have 
been compared to each other and to the fields given in the present FCC 
Ground Wave Field Intensity curves. 


SURVEY OF METHODS FOR PREDICTING GROUND WAVE FIELD INTENSITIES 


A bibliography is given at the end of this report listing publications con- 
sulted in the preparation of this report. Bracketed numbers in the report 
refer to the same numbered reference in the bibliography. It appears that 
there is reasonably close agreement in the reference publications for 
prediction of ground wave field intensities for distances of about 100 km 
or less. For this distance range predictions are usually made using the 
Sommerfeld (1) surface wave field equations which were also treated by 
Norton (4), (5), and (7). Surface wave fields are predicted by use of 
the complex error function as discussed by Norton. At greater distances 
diffraction must be taken into account, and different calculation methods 
are required. For these greater distances, predictions are usually made 
using a residue series originally developed by Watson (2). Due to slow 
convergence of this residue series at smaller distances, there is a range 
of distances where neither the surface-wave equations nor the residue 
series, unless carried to a great number of terms, accurately predict 
ground wave field intensities. 


Various authors have used different expedients to bridge this gap. Norton (5) 
suggests a method earlier proposed by Eckersley (6) to draw a continuous field 
intensity curve versus distance. The method consists in making two graphs of 
field intensity versus distance using identical scales of distances on the 
abcissae and field on the ordinate. One graph contains calculated surface- 
wave field intensities; the second, calculated residue-series field inten- 
sities. The graphs are superimposed on each other on a light table and the 
distance scales aligned. The graphs are moved vertically until the two curves 
are tangent to each other. 


A composite curve is then prepared by using the surface-wave fields to the 
point of tangency and adjusted residue series fields beyond this point. The 
residue series adjustment consists in multiplying these field values by the 
ratio of the surface-wave field to the residue field at the tangent point. 


Another method of preparing a composite ground-wave field-intensity curve 

has also been suggested by Norton (7). In this second method, Norton suggests 
that the surface-wave equations be used up to distances of 50/ FMhz 1/3 miles 
to prepare this part of the field intensity curve. At a considerable distance 
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beyond the optical horizon а second curve is constructed using equations 
and graphs given in this reference. The graphs and equations relate field 
intensities to ground constants, frequency, and the radius of the earth 
as modified by diffraction. Between these two curves a smooth curve is 
drawn. 


Bremmer (13) has suggested a third method of deriving a ground-wave field- 
intensity curve. He suggests that the Norton-Sonmerfeld surface-wave field 
equations be modified to include а spherical-earth correction factor. (The 
Sommerfeld surface wave was derived for plane earth). The correction factor 
is a function of the effective earth radius, the wavelength, and earth 
constants. Bremmer states that with this correction the surface-wave fields 
will merge with the residue series fields approximately at the effective 
optical horizon. This merger may possibly occur for frequencies considerably 
above those for the Standard Broadcast band. However, tests of the Bremmer 
equations for Standard Broadcast frequencies with usual values of ground 
constants show that the residue fields are almost always smaller than the 
surface fields and do not merge with either the corrected, or uncorrected, 
surface-wave fields. 


In addition to the fact that the surface-wave and residue-series curves do 
not merge, an even more serious problem, which has been encountered with the 
Bremmer spherical-earth correction to the surface-wave fields, is that the 
corrected surface-wave fields are erroneously large for high conductivities 
such as sea water at Standard Broadcast frequencies. 


Bremmer does give considerable information on the residue series in this 
publication, and it is possible, if desired, to derive an unlimited number 

of residue-series terms from the equations given by Bremer. Bremmer states 
that only three terms of the residue series are needed for most calculations. 
Al'Pert (11) used only two terms of this series in his field predictions. 
However, computer test programs written as part of this study show that at 
least six terms of the series are needed for convergence to 1 percent of the 
ultimate field value for distance of 100 kilometers at broadcast frequencies 
for the range of ground constants which may be encountered. More terms are 
needed at smaller distances. It has been found during the course of computer 
testing the Bremmer equations, that the uncorrected surface-wave field graphs 
parallel the residue-series field graphs for distances ranging from арргохіта- 
tely from 40 to 110 kilometers. This parallelism makes a smooth merging of 
these curves possible by readjusting the ordinate values of the residue fields 
as suggested by Norton (5) with regard to merging of the surface-wave fields 
and Watson fields. If the ordinates of the residue series are not adjusted in 
this manner, a reverse "S" joining curve will be necessary to connect the two 
parallel primary curves. Both methods of joining the two fields have been 
tested. It appears that the method of readjusting the residue-series ordinates 
gives a smoother curve. 


ата 


COMPUTER PROGRAMS DEVELOPED FOR GROUND WAVE FIELD INTENSITY PREDICTION 


Three different computer programs have been developed to predict ground-wave 
field intensities using the three methods discussed in the previous section. 
The programs are built up from a number of subroutines chosen to perform 
the various calculations required. This type of program construction makes 
it easy to modify programs by changing only a particular subroutine without 
altering other parts of the program. Since all three programs calculate 
fields within the optical horizon using the Norton-Sommerfeld equations, 
only one set of subroutines is used by all three programs to calculate 
fields within the horizon distance. Beyond the optical horizon, the three 
programs differ in computation methods. Two of the programs, called гезрес- 
tively the Watson-Sonmerfeld and the Bremer, use different versions of the 
Watson residue series to calculate fields. The other program called the 
Norton 1941 uses a different set of formulas for this purpose. 


Ultimately only one ground wave field prediction program is needed. However, 
during program development it has been found useful to have three different 
programs for cross—checking purposes. Fields predicted by these programs have 
also been compared against fields given in the FCC Ground Wave Propagation 
Curves. Figures 1-4 compare predictions by these four sources at .55 and 
1.6 MHz for conductivities of 2 millimhos per meter and 5000 millimhos per 
meter and relative dielectric constants of 15 and 80 respectively. Except 
for Figure 3, these figures show very close agreement between the fields 
taken from the FCC Ground Wave Field Intensity Curves and those predicted 

by the Norton-Sonmmerfeld surface-wave computer program at distances within 
the effective horizon distance. Figure 3 shows the Norton-Sonmerfeld fields 
falling several percent below the FCC curve values for distances within the 
horizon. It is believed that this difference is caused by the fact that the 
Norton-Sommerfeld fields in this figure were calculated for a frequency of 
1600 kHz whereas the FCC fields which were taken from the FCC 1520-1600 kHz 
curves were probably calculated for a frequency of 1560 kHz, the mid fre- 
quency fer this set. 


Comparison of the fields shown for a conductivity of 2 millimhos, within 
horizon distances as given by the FCC 1430-1510 kHz curve and the FCC 1520-1600 
kHz curve, show that these fields are approximately proportional to the inverse 
frequency ratio squared. On this basis, it is to be expected on Figure 3, that 
the Norton-Sommerfeld calculated fields would be about 5 percent less than the 
FCC fields within the horizon. This difference between the calculated fields 
and the FCC fields does not appear on Figure 4 because for sea water both fields 
are so near to the inverse distance field that there is very little effect from 
a variation of frequency. 


The close agreement for within horizon distances between the FCC fields 
and the Norton-Sommerfeld fields shown in Figures 1 and 2 probably results 
from the fact that both curves were calculated for the same frequency of 
550 kHz. The trends shown in these four figures have also been found to 
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apply generally in numerous other comparisons. That is, for distances 
within the effective horizon, there is good agreement between FCC field 
predictions and those generated by the computer programs. The agreement 
among these four sets of fields is poorer over the horizon, and none of 
the computer programs gives fields which consistently agree with the FCC 
fields. 


DESCRIPTION OF COMPUTER PROGRAM LISTING 


Included in this report is a listing of the Fortran programs used during 
this study for calculation of ground wave field intensities. Although a 
MAIN program is included in the listing, various other versions of the 
MAIN program are possible depending upon which sets of calculations are 
desired. The function of the MAIN program is to merely introduce the 
variables—--frequency, ground constants, and distances-—-to the various 
subroutines and then write out the calculated results in some desired 
format. 


No large effort has been expended in the program development to choose the 
most efficient calculation methods, and undoubtly, some reduction in calcu- 
lation times can be made by rearrangement of some of the subroutines. Рог 
example, many parameters which are functions of frequency and ground соп- 
stants alone, but not of distance are calculated in the subroutines. It 
is therefore efficient to calculate these parameters once for a large set 
of distances. BREMR and SBREMR illustrate this arrangement by only calcu- 
lating the residue series coefficients when frequency or ground constants 
are changed. A similar bypass for redundant calculations could probably 
be incorporated in subroutine ABG. 


The three methods of calculating ground-wave field intensities used in this 
study as stated previously are called Norton 1941, Bremmer, and Набзоп- 
Sommerfeld respectively to indicate the source of each program. All three 
programs calculate fields within the optical horizon using the same set of 
subroutines КОМО and FPB. The three differ in calculating over-the-horizon 
fields. NUMD calculates P, the numerical distance, and B, the phase angle 
of P, from the ground constants, the frequency, and the distance. 


Subroutine FPB then uses P and B from NUMD to calculate A, the surface wave 
attenuation factor below inverse distance fields, and PHI, the phase angle 

of A. FPB uses three different methods to calculate A depending upon the 
value of Р. The methods were chosen to reduce calculation times by тіпіті?- 
ing unnecessary iterations. For values of P equal to, or less than “7%, 2 

is calculated by use of the М function (14, р. 297). For values of P between 
.Т and 5, FPB uses a convergent infinite Р series given by Norton (7, р. 637). 
For P values exceeding 5, A is calculated in a closed form using an equation 
for the W series (14, p. 328). 
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Taking the three methods in order, the Norton 1941 program calculates оуег- 
the-horizon fields using the subroutines FMVMA and ABG. FMVMA constructs a 
transition curve beyond the fields within the optical horizon and the fields 
at points well over the horizon using NUMD, FPB, and ABG. ABG calculates 
fields over the horizon using equations given by Norton in reference (7). 


The second program called Bremmer uses residue series equations given by 
Bremmer (13) to calculate over—the-horizon fields. It uses subroutines 
FMVMB, BREMR, and SBREMR for these calculations. It has been found by 
experiment that the Bremmer residue fields lie in curves that parallel 
the surface wave fields curves for an appreciable distance on either side 
of the horizons. ЕМУМВ and the MAIN program use this parallelism to join 
these curves at the horizon by multiplying the residue series fields by 
the ratio of the surface wave field to the residue field at the horizon. 


The third program called the Watson-Sommerfeld program is based on formulas 
and tabular data given by Norton 


The Watson-Sommerfeld program calculates over—the-horizon fields using the 
Watson diffraction formula corrected for refraction in the lower atmosphere 
as suggested by Burrows (9), and incorporating a correction by Eckersley (6) 
for finite earth conductivity. Fields given by the Watson diffraction Гог- 
mula are less than the surface wave fields at equal distances. Consequently 
to join these two fields into a continuous curve it is necessary to move the 
Watson curve vertically until is is tangent to the surface wave curve. 
Watson fields are calculated by subroutines FOH and SFOH based on Norton (4). 
The joining of the two curves, which was formerly accomplished mechanically 
be superimposing graphs of surface-wave fields and Watson fields and then 
moving the Watson graph vertically for tangency, is now performed by sub- 
routine MERGE. MERGE finds the distance at which the slopes of the two 
field curves are most nearly equal and sets a multiplier of the Watson 
field to produce equality of the two fields at this distance. Fields are 
calculated for distances beyond the tangent distance using the adjusted 
Watson equations. The surface wave equations are used up to the tangent 
distance. 


Equations and tabular data on which the various subroutines are based are 
given in detail in the various references and іп the interest ОҒ brevity 
are not repeated here, since the equations and tables are also given in 
the Fortran listing of these programs which is included in this report. 


абы 


DISCUSSION AND RECOMMENDATIONS 


During this study numerous comparisons have been made of the fields рге- 
dicted by the listed computer programs and the fields predicted by the 
FCC Ground Wave Field Intensity curves. These comparisons have shown 
the same general trends evident in Figures 1-4 of this report, that is, 
there is good agreement between the computer program field predictions 
and the FCC curve fields out to the optical horizon at distances between 
80 and 100 kilometers. At greater distances, field predictions of the 
computer programs and the FCC curves differ appreciably. A major 
influence on this difference appears to be the ground conductivity. 

For example, for sea water the Norton 1941 program predicts fields which 
are nearer to FCC curve values than those of the other two programs. For 
this conductivity the Norton 1941 program predictions are usually within 
25 percent, or less, of the FCC curve fields. For lower ground сопдис- 
tivities, the field predictions of the other two programs more nearly 
approach the FCC curve fields than does the Norton 1941 predicted fields. 
The Bremmer program usually predicts fields which аге about 10 percent 
above the FCC fields for distances of 100 to 200 kilometers, and then 

at greater distances falls more and more below the Fcc fields to a value 
about 50 percent lower than the FCC fields at 1000 kilometers. The 
Watson-Sommerfeld program fields have an opposite slope to that of the 
Bremer fields having roughly the same positive error difference in the 
100 to 200 kilometer range and then the Watson-Sommerfeld fields continue 
to increase with distance above the FCC fields to a value about 50 percent 
above the FCC fields at 1000 kilometers. 


It should be noted in passing that, which, if any, of these four sources 
of fields predictions is correct is not known. The region where the major 
differences of predictions lies, from 100 to 1000 kilometers, is a region 
where very few sets of field measurements have been taken. There are many 
reasons for the lack of measured ground wave field data in this region. 
Among these reasons are the following: (1) Field measurements for most 
broadcast station proofs of performance are made at distances which are 
less than 100 kilometers, (2) Except for all-sea-water paths, most other 
paths for distance of 100 kilometers, or more, will usually cross several 
regions of different conductivities making it difficult to relate ground 
constants and fields, (3) Only around mid-day are skywave signals not 
present over most of this range complicating determination of the true 
ground wave field values, (4) Interference from co-channel and adjacent 
signals may contaminate the desired ground wave signal, (5) At the larger 
distances man-made and atmospheric noise may have magnitudes which are 
comparable to the ground wave signal. 


гд 


In view of the previous discussion which shows that, first, there is not 
consistent agreement between any of the computer predicted fields and the 
FCC curve fields, and second, that it is unlikely that field measurements 
сап be made to find which method of predictions is correct, there appears 
to be no obvious choice of a method which is consistently better than the 
other methods for revising the FCC Ground Wave Field Intensity Curves. It 
would appear that complications which may arise due to revision of these 
curves should be weighed against advantages which may be obtained from 
having a computer program which gives consistent and reproducible field 
values at all distance and ground conductivities. 


Of the computer programs developed in this study it is believed that the 
Bremmer program has the best theoretical justification being based on the 
more recent theoretical work and providing sufficient terms of the residue 
series to provide any degree of convergence desired in over-the-horizon 
ground wave fields. Within the effective horizon there is no difference 
between the programs since they use the same subroutines for field calcu- 
lations. 


Accordingly, if it is desired to choose a computer program from the three 
listed programs for revision of the FCC Ground Wave Field Intensity Curves, 
the Bremmer program is recommended. It is further recommended that some 
additional effort be expended in finding methods to increase the саїсціа- 
tion efficiency of the program by eliminating all redundant operations. 
Calculation of the residue series fields and the surface wave fields 
involves mathematical operations with complex numbers. It is possible 
that some increase in program calculation efficiency would result if 

the programs were changed to use the 1106 computer complex number routines 
rather than those presently incorporated in the subroutines. It is also 
suggested that the Bremmer program be used to prepare tabular arrays to 

be used іп а table-look program to be developed which will have a faster 
execution time than the Bremmer program. Not considered in this report 
are two ground-wave field intensity programs in development by ITS and 
CCIR, since they were unavailable for use during this program. However, 
an interim version of the CCIR program was studied in this project. This 
program also gave fields at variance with the FCC curves. It is understood 
that final versions of these two programs may be available in 1979. 


The present propagation programs do not provide several procedures which are 
often needed in ground wave propagation studies. Among these missing ргосе- 
dures are: (1) A capability, given a value of field intensity, to find the 
corresponding distance, and (2) A capability to find field intensities for 
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ground маме paths with multiple conductivities. Neither of these procedures 
is considered difficult to implement with conventional computer programming. 
The first procedure probably could be implemented most efficiently if the 
present propagation programs were used to prepare tables of field intensity 
versus distance for the required values of ground constants and frequency. 
The distance for a given field intensity would then be found with a table- 
look-up subroutine which would enter the tables on the distance axis and 
look for the given field intensity. The subroutine would interpolate 
between tabular fields when the desired field was not a tabular value. 

This first capability could also be provided, somewhat less efficiently, 
with a subroutine which would use the present propagation programs. This 
subroutine would assume distances in an iteratively, convergent manner, 
checking the desired field and the field at the assumed distance and then 
modifying the distance so as to bring the two fields nearer together until 
the desired degree of convergence was obtained. 


The second capability, to predict ground wave fields for paths with multiple 
conductivities, can easily be programmed using the previously discussed 915- 
tance/field subroutine and the existing propagation programs. The subroutine 
would use the present propagation program to find field at the first conduc- 
tivity boundary. The subroutine would then use this field with the distance/ 
field subroutine to find the equivalent distance for this field for the 
second conductivity. The subroutine would ада the path length for the 
second conductivity to the equivalent distance and use the present propaga- 
tion program to find the field at the next conductivity boundary. The sub- 
routine would repeat this process as many times as required for the various 
conductivities. 
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